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Abstract 



Spectral anomaly for interacting Fermions is characterized by the spectral 
function A([k — Icf], w) satisfying the scaling relation A(A Vl [k — /cp], A^u;) = 
A yA A([k — kp],uj), where yi, 1)2, and z/a are the exponents defining the uni- 
versality class. For a Fermi liquid y\ = 1, j/2 = 1, DA = — 1; all other values 
of the exponents are termed anomalous. In this paper, an example for which 
Vi = 1) 2/2 = 1) but yA = 01 — 1 is considered in detail. Attractive interaction 
added to such a critical system leads to a novel superconducting state, which 
is explored and its relevance to high temperature cuprate superconductors is 
discussed. 
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I. INTRODUCTION 



The normal state of high temperature superconductors is unusual: although the normal 
state has a metallic response and reflects gaplessness, its properties cannot be explained by 
Fermi liquid theory However, a microscopic theory of such a gapless non- Fermi liquid 

in dimensions d > 2 does not exist. In contrast, the breakdown of Fermi liquid theory in 
d = 1 has been known since the pioneering work of Mattis and Lieb [ffl. The absence of 
fermionic quasiparticles, spin-charge separation, and power-law correlation functions have 
been extensively discussed Pfl, and Haldane J7[ has shown that these properties hold for a 
large class of interacting systems, the Luttinger liquids. 

At the same time, a theory of high temperature superconductors has been proposed |||9j 
in which Josephson tunneling of the Cooper pairs between the CuO layers plays an important 
role, and leads, in particular, to a striking dependence of the transition temperature on the 
number of such layers in the unit cell and to a high transition temperature by exploiting, in 
a very special way, the new energy scale proportional to the square of the tunneling matrix 
element divided by the in-plane band width. The validity of this theory depends on the non- 
Fermi liquid nature of the normal state |H| , as the special momentum conserving Josephson 
pair tunneling Hamiltonian cannot possibly be correct for a Fermi liquid. In view of the 
absence of a precise microscopic theory of the non- Fermi liquid, it is necessary to explore the 
extent to which a consistent theory can be developed based on general arguments because 
the consequences of this interlayer tunneling theory are quite tantalizing. It is our intention 
to present some general arguments in a reasonably complete form and to point out directions 
in which further improvements are essential. 

In the present paper, we shall assume the existence of a higher dimensional non-Fermi 
liquid state and study the consequences for superconductivity. The essence of our arguments 
is a model of the one-particle Green's function and the corresponding spectral function. In 
constructing this model, we follow the stratagy of Galitskii and Migdal [|13[ who constructed 
a model for a Fermi liquid and examined its self consistency. The resulting theory is, of 
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course, the same as the Landau theory of Fermi liquids. In the Galitski-Migdal model, the 
one-particle Green's function, analytically continued to the complex plane, contains simple 
poles on a second Riemann sheet, while the principal sheet is cut along the real frequency 
axis. The quasiparticle is a resonance in a many-body state, resembling the bare particle. 
In contrast, in our model, there are no poles but branch points on a second sheet, located 
infinitesimally close to the real axis. Thus, the non-Fermi liquid state is gapless and responds 
like a metal, but without quasiparticles. Such a model is motivated by the Luttinger liquid 
concept, but the spectral function is only qualitatively similar to that of the Luttinger 



liquid. In fact, in the present context, this model was briefly sketched by Wen pl| . The 
key element is the characterization of the non-Fermi liquid state as a critical state defined 
by a scale invariant spectral function ||10|| . This characterization is described by the term 
spectral anomaly, and, as in the theory of critical phenomena, the universality classes are 
distinguished by the critical exponents of the spectral function. Only a very special choice 
of the exponents leads to the fermi liquid model. In fact, we believe that all systems that 
respond as metals can be be described by our model. We also pay special attention to 
the Ward identity that relates the vertex function to the inverse Green's function. This is 
essential in building a consistent description of the non-Fermi liquid state. 

The outline of our paper is as follows. In Sec. (II), we introduce the concept of spectral 
anomaly and examine the nature of the normal state. In Sec. (Ill), we consider the effect 
of various interactions added to the normal state described by spectral anomaly. In this 
section, we discuss only in-plane pairing interactions. In Sec. (IV), the superconducting gap 
equations with in-plane pairing interactions are solved for both s-wave and d-wave kernels, to 
examine the effect of the non-Fermi liquid behavior of the normal state. Sec. (V) contains a 
discussion of the Josephson effect and consititutes a qualitative justification of the interlayer 
pair-tunneling Hamiltonian discussed in Sec. (VI). It is worth noting that the analysis in 
this section bears some resemblance with a recent paper by Khlebnikov jnj. In Sec. (VI), 
we discuss the interlayer tunneling model, and, compute tunneling density of states with 
parameters as realistic as possible to compare with recent vacuum tunneling experiments. 



In Sec. (VII) we consider spin-charge separation and discuss how our basic model is changed 
in its presence. In the final concluding section, Sec. (VIII), we point out the directions in 
which our theory should be extended and improved. 

II. SPECTRAL ANOMALY 

For a finite system, the Kallen-Lehmann representation fl4| of the time-ordered one- 
particle Green's function exhibits a discrete set of poles corresponding to the excitation 
frequencies of the system. In the infinite volume limit, the poles merge into a cut that runs 
along the real axis. Moreover, causality requires that the analytic continuation from the 
real axis to the complex plane be free of singularities in the first and the third quadrants. 
In fact, the only singularity on the principal sheet is a cut along the real axis. 

To describe a Fermi liquid, Galitskii and Migdal |13[ proposed a model of the one-particle 
Green's function, which contains singularities in the second and the fourth quadrants on a 
second sheet. These singularities were chosen by them to be simple poles, sufficiently close 
to the real axis to define long-lived quasiparticle states. The pole in the fourth quadrant 
corresponds to a quasiparticle and that in the second quadrant to a quasihole. (Throughout 
the paper we choose to work with a fixed chemical potential fi in a grand canonical ensemble 
rather than in a canonical ensemble for formal simplicity). This model corresponds to the 
assumption that the matrix elements of the Fermion operators between the ground state and 
the exact eigenstates are sharply peaked at an energy that depends on the wavevector k. 
The quasiparticle pole in the complex plane is the representation of the physical phenomena 
on the real frequency the theory of analytic continuation implies. 

A simple pole corresponds to a unique energy-momentum dispersion of the quasiparticle: 
with a wave vector k, we associate an energy e^. And, with the renormalized wave function, 
the quasiparticles can be made to resemble the bare particles for energies asymptotically close 
to the Fermi energy. Other properties of the quasiparticles also follow from the assumptions 
of the model. The charge of the quasiparticle is unity and its spin 1/2. 
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The Galitskii-Migdal model can be summarized by the spectral function 

A(k,u) = --ImG R (k,uj) 

71 

= z k 8{u-{e{k)-n)), (2.1) 

where Gr is the retarded Green's function, and z k is the quasiparticle residue. As this model 
is valid for quasiparticles close to the Fermi surface, we can linearize the spectrum and write 
e(k) — /i ~ vp{k — hp), where kp is the Fermi wave vector and vp is the Fermi velocity. 
Then, the spectral function has the scaling property 

A(A[k - k F ],Auj) = A^A^k - k F ],u), (2.2) 

reflecting the gaplessness of the elementary excitations. The system is poised at criticality. 
That a normal Fermi system represents a critical system is, of course, understood. The im- 
portant point, however, is that the criticality in a normal Fermi system does not usually lead 
to serious dynamical consequences. This is because the modes (the Fermionic excitations) 
decouple, as in the Gaussian model of classical critical phenomenon. 

We would like to build a model in which the criticality of the Fermi system is more 
important. The simplest possible hypothesis is to assume that the analytic continuation of 
the Green's function to a second sheet contains branch points instead of simple poles. We 
postulate, therefore, a spectral function that satisfies the scaling relation 

A(A vl [k-k F ],A V2 u)=A VA A([k-k F ],u), (2.3) 

where y±, y 2 , and yA are the exponents defining the universality class of the critical Fermi 
system. The values of the exponents other than the set y± — 1, y 2 — 1, and yA = — 1 
will be termed anomalous, while this special set represents a Fermi liquid for which the 
branch points collapse into simple poles. The notion of a spectral function with anomalous 
exponents will be termed spectral anomaly. 

The anomalous spectral function implies that a given k does not correspond to a single 
frequency but a continuum of frequencies. The creation operator of wave vector k, applied 
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to the ground state, creates a state that couples essentially to many eigenstates that cannot 
be expressed by the model of a weakly damped quasiparticle. 

To establish continuity with the Fermi liquid model, as the exponents approach the 
Fermi liquid values, we make an important assumption. The location of the branch point 
is assumed to be in one-to-one correspondence with the states of the noninteracting system, 
which may be thought of as a remnant of Luttinger's theorem. This is by no means an 
obvious assumption, but holds for one-dimensional electron gas models [§J, and also for 
gauge models [|l5j in higher dimensions. 

From above, it is clear that we have not incorporated spin-charge separation in our 
definition, for we have assumed that the spectral function is singular at a single energy. The 
mathematical signature of spin-charge separation can be understood from the properties of a 
one-dimensional Luttinger liquid where the electron is a composite particle. The elementary 
particles are spinons that carry only spin and ho Ions that carry only charge. As a result, 
the spectral function is singular at two energies. A generalization to include spin-charge 
separation is discussed in Sec. VII. 

From the dispersion relation we can now determine the real part of the retarded Green's 
function, i.e, 



ReG R {[k-k F },u) = / 



P r+oo 



dw 



,lmG R ([k - k F ],u) 



(2.4) 



uj — to' 



Thus, the real part also satisfies the same scaling relation given by 



ReG R (A yi [k - k F },A y2 u) = A VA ReG R ([k - k F ],u). 



(2.5) 



From the generalized homogenity assumption, we also get 



A{[k-k F ),0) 



1 



A([k - k F ]/\k - k F \,Q) 



(2.6) 



k - k F \ VA f yi 



and 




(2.7) 



S 



It is also possible to derive the momentum disrtibution function without any detailed knowl- 
edge of the function A(k, uj). From the Lehmann representation, the momentum distribution 
function n(k) is given by 

f o 

n(k-k F ) = duA([k-k F ],u). (2.8) 

J — oo 

Using the scaling relation, we can write 

n(k - k F ) = ( duoA~ yA A{A yi [k - k F ], A V2 uj). (2.9) 

J —oo 

Because A is an arbitrary scale factor, we can chose A yi \k — k F \ = 1. Then, 

n(k — k F ) — \k — k F \ {y2+yA)/yi f dxA(Sign(k - k F ),x). (2.10) 

J — oo 

The discontinuity present in a Fermi lquid is destroyed due to spectral anomaly. In a Fermi 
system, the occupation of a state, n(k), cannot diverge, and the inequality (t/2 + yA)/yi > 
must be satisfied. There is a superficial difficulty with the scaling argument above: the 
integral for n(k) does not converge. Thus, the quantity n(k) cannot be obtained from the 
scaling argument, and we need to restore the cutoff. However, we believe that the critical 
exponent is still given correctly by the scaling argument. 

Consider now the time ordered Green's function, G(k,cu). It is 

G(k, u) = ReG R (k, uj) - iSign(co)lmG R (k, to). (2.11) 

In terms of the time-ordered Greens function, we can write 

— j r+oo 

n (k) = — lim / dwG(k,uj)e luJt . (2.12) 



This integral can be broken up into three pieces. We write 



n(k) = — lim 

V ^ 27Tt-0+ 



/— u) c r+Uc r+oo 

duG{k,u)e twt + / duG(k,u)e tuJt + / duG{k,uj) 
-oo J— U) c J +w c 



e iwt 



(2.13) 



The Green's function is assumed to be anomalous in the frequency range between — uo c and 
uj c . In the second integral, we can set t = and apply the scaling form. The integral is well 
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behaved becuse the singularity, as uj — > 0, is softer than the Fermi liquid case (oc ^). The 
first and the third integrals can be estimated by replacing G by its correct asymptotic form, 
which is G(k,cu) — > ^, as uj — > oo, resulting from the Fermion anticommutation relation. 
Finally, we find again that 

n{k) = \-C Sign( k -k F )\k-k F \ ? ( 2 . 14) 

where C is a constant. 



A. The basic model 

To proceed further, we introduce an explicit model that is defined by the retarded Green's 
function: 

Gr(*, u>) = -a; c < w < Wc . (2.15) 

where is the phase, a is the anomaly exponent, and 5 is an infinitesimal positive number. 

To conform to causality, the upper half plane should be free of singularities, and the cut 

is assumed to lie in the lower half plane. The quantity g(a) will be determined from the 

normalization. In a rough manner, one may include all nonsingular self energy effects in the 

definition of e(k). Therefore, 

1 sin0 6(w-e(k)) sin[0 - tt(1 - a)} 9(e(k) - u) 
ImGR(/c,cj) = — g(a) -. — ■ —. g(a) -ttt\ — . 

7T ^ 7T U)«(U0 - E{k)y- a ' 7T UJ^(e(k) - Uj) l ~ a 

(2.16) 

Because the spectral function has to be positive, the phase must satisfy the inequality 
> > — Tra. This condition is not restrictive enough to fix uniquely. However, the time 
reversal symmetry is violated unless = —an/2. This is because of the invariance with 
respect to uj — > — uj, and e(k) — > —e(k). Therefore, 

A(k,uj)=g(a) y2 \ — (2.17) 

v 1 y 1 tiuj? \e(k) -uj !- a v y 
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When e{k) = 0, our spectral function is similar to that of a Luttinger liquid without spin- 
charge separation in d — 1 |I6| , |T7| , but it is very different when e(k) ^ 0. In a Luttinger 
liquid, the spectral function diverges with a power law as uj — > vpk from above, vanishes 
for —vpk < uj < vpk, and vanishes with a power law as to — > —Vpk from below. In fact, 
the similarity with the one-dimensional system cannot be expected. The special kinematic 
constraints present in d = 1 are missing in higher dimensions. Finally, note also the difference 



with the gauge models [ 15|| : although the spectral functions are similar for e(k) = 0, for 
k ^ 0, the self energy in the gauge models is nonsingular. 

To preserve the equal-time anticommutation relation of the Fermions, one must satisfy 
the sum rule Jf^ dujA(k, uj) = 1. If we integrate the spectral function, given in Eq. ( p,17|) , 
we get instead: 

<k)Y . A e(k)Y 



sinf 

duM(fc,u;) = ^(a) — ^ 
7ra 



1 + 



(2.18) 



If we choose 



7ra 



»(«) = 9 • ( 2 - 19 ) 
2sni(y) 

the anticommutation relation will be satisfied if \e(k)\ u c , which is precisely the regime 
in which a scaling theory is appropriate. 

We are still left with a minor problem with respect to the choice of the phase. For the 
exact retarded Green's function, the phase at uj = oo, 0(oo) = — e, because the imaginary 
part of Gr vanishes faster than the real part. As uj runs along the real axis, IhiGr does not 
change sign and the phase remains in the lower half plane, reaching </>(— oo) = tc + e. For 
our Green's function, the phase remains in the lower half plane, varying between —an/2 
and (tt + air/ 2). Clearly, the exact limits at ±oo are violated, unless the scaling form is 
supplemented by the exact asymptotic limits of the Green's function. However, as long as 
a 7^ 0, discontinuities of the phase remain in our model; however, these disconitinuities 
appear to be physically inconsequential. 

The time-ordered Green's function, G(k,uj), is given by 
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G{k,u) = g(a)uZ 

Note that the argument of the ^-function is the frequency variable, u, and not e(k). For a 
Fermi liquid, there is no difference between these two cases; this is not so here. The correct 
argument is always uj (the total energy), and not the kinetic energy e{k). 



1 2 



9(u) 



(u - e(k) + i6y- 



+ 



1 2 



[uj - e(k) - iS) 1 - 



[2.20) 



B. Momentum distribution function 

The momentum distribution function is now given by 

a f° 1 
n(k) = --— — / duo— — . (2.21) 

We get 

(2.22) 

To understand what we meant previously by cutoff effects, let us attempt to derive n(k) by 
scaling and extending the limit to — oo. Then, 

n(k) = | 

We see that the singularity of n(k) is captured correctly, but not the correct functional form. 
This is as it should be because the remaining integral does not converge. The calculation of 
n(k) from the time-ordered Green's function gives, as before, the same result. 



n{k) 



1 

2 



1- Sign (£(£)) 



e(k) 



e(k) 



0J r 



dx-. — 

oo \X 



Sign(e(fc))| 



l-a 



(2.23) 



C. Matsubara Green's function 

The Matsubara Green's function, G(k,iv n ), can be obtained from the spectral function. 
Note that 

Q{k,u n )= r du'^1, (2.24) 

J-oo tUJ n — UJ' 

where u) n = n{2n + 1)T are the Matsubara frequencies; T is the temperature. Substituting 
the spectral function, we get 
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Q(k,uj n ) = g(a 



where we have used the integral 



9(u n )e l ™ + 9(-u n )e 
^{iu n - e{k)y- a 



2 



^2.25) 



r.a-1 



Ti 



f°° d X " 
Jo (z + x) smina) 



_a-l 



(2.26) 



with |Arg(z)| < 7r and 1 > Rect > 0. 

D. Density of states 

The density of states, n(u), can be calculated from the spectral function. It is 



n(uj) 



a r dJ p{u>) 



2ui 



\uj' — UJ 



l-a' 



(2.27) 



where 



k 

If we choose p{uS) ~ p = const, for — W < u <W and otherwise, then 



^2.28) 



W v " " 



W _ LU\ a ' 
0J C UJ C 



(2.29) 



We have also numerically calculated n(u>) for a more realistic p{uj). For small a, the smearing 
is not enough to destroy the van Hove singularities, while for larger a it can completely 
smooth out the singularities. 

In contrast to Luttinger liquids, the density of states does not vanish as uj — > 0. In 
fact, it is essentially unchanged from the Fermi liquid value, obtained as a — > 0. It would 
agree exactly with the Fermi liquid value if uj c is chosen to be W . Once again, the correct 
treatment of the cutoff cannot be overemphasized. The scaling trick, applied to the integral 



in Eq. ( ^27|) leads to 



n(u) 



pa 


UJ 


a roo 


T 




J — oo 



dx- 



x — Sign(cu) 



l-a' 



^2.30) 
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where the limits of the integral were set to ±00. It would appear that the density of states 
vanishes as uj —>■ 0. This is incorrect, because the remaining integral does not converge; the 
cutoffs cannot be set to ±00. The vanishing of the density of states at the Fermi energy in 
a one-dimensional Luttinger liquid is a result of the special kinematic constraints that do 
not easily generalize to higher dimensions. 



E. Consistency of the basic model 

It is useful to consider Fermi liquid theory as an effective low energy theory of metals 



IgJlq-pOf . We shall briefly repeat the arguments due to Polchinski | 20|j . Imagine that we 
are considering an effective field theory below an energy scale uj c at which there are strong 
Coulomb interactions. We want to find the effective low energy theory for energies smaller 
than uj c . The elementary excitations of the effective field theory underlying the Galitski- 
Migdal model are charged, spin-1/2, quasiparticles, similar to the original electrons. 
Consider first the free action: 

J dtd d p{ir p At)dtM) - (e(p)-fi)r p At)M)} (2-31) 

As we scale all energies by a factor s < 1, the momenta must scale to the Fermi surface. To 
accomplish this scaling, we write 

p = k + l, (2.32) 

where k is a vector on the Fermi surface and 1 is a vector orthogonal to it. Let 

e(p)-fi = lv F (k) + 0(l 2 ), (2.33) 

where t>F(k) is the Fermi velocity. As we scale the energies and the wave vectors, E — > sE, 
k — > k, and 1 — > si, the action remains fixed, provided the dimension of the fermion operators 
is s -1 / 2 . Consider now all possible terms consistent with the symmetries and examine their 
scaling. If there are no relevant operators, the effective low energy theory is self-consistent. 
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(The generated mass term can be absorbed by redefining the fermi surface to be the true 
interacting Fermi surface) 

The most important interaction is the four fermion operator: 

J dtd^ 1 ^ dl l( l d - 2 k 2 dhd^ks dl^" 1 ^ dhV{k u k 2 , k 3 , k 4 ) (2.34) 

r P1 At)^ P3 At)r P2 At)^(t)5 d ( Pl + P2 - p 3 - p 4 ), 

where we have assumed that the interaction is short-ranged. It can be seen that the interac- 
tion scales as s times the dimension of the delta function. Let us assume, for the moment, 
that 

5 d (pi + p 2 - Ps - P 4 ) = 5 d (ki + k 2 - k 3 - k 4 + li + h - 1 3 - U) 

-^(kx + ka-kg-k^, (2.35) 

where we have ignored 1 in comparison to k. This justified because the vectors / scale to 
zero. Now the argument of the 5-function does not depend on s, and the four fermion 
interaction scales as s, vanishing in the limit s —>■ 0. The interaction is irrelevant. The 
conclusion is almost correct except for the special kinematics implied by the Cooper channel 
scattering processes. For these processes, the delta function can be seen to be of order 
transforming the dimension of the interaction to s . The net result is that the interaction is 
marginal pOfl . Thus, except for special kinematics, the four fermion interaction is irrelevant. 
The argument, valid in dimensions d > 2, fails in d — 1. In d — 1, the 5-function is always 
of order s -1 and the interaction is always marginal. 

We now turn to the consistency of our basic model. The action is 

J dud d pG-\p,u;)r p M^ P ^), (2-36) 

where the Green's function G corresponds to that of the basic model. The dimension of 
V> CT (p, uj) is s~( 3_a )/ 2 and that of ^> CT (p, t) is s _ ( 1_Q )/ 2 . If we follow Polchinski, we see that the 
the four fermion interaction is irrelevant; even for those exceptional kinematics for which the 
interaction was marginal in the Fermi liquid case, it now scales as s 2a . The spectral anomaly 
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is more stable than the Fermi liquid! In fact, in the weak coupling regime, it does not even 
allow a superconducting instability. In Sec. V, we shall see that this is indeed correct. The 
coupling has to reach a threshold before superconducting instability occurs. Remarkably 
enough, the Josephson pair tunneling Hamiltonian || describing interlayer tunneling is 
relevant and results in a superconducting instability for arbitrarily weak coupling. Such 
a term must therefore be included in the effective low energy theory. 



F. Vertex function 



In the normal state, the charge vertex and the current vertex satisfy the Ward identity 



21 



£ q^(k + |, k - f ) = G-\k + |) - G~\k - f ). (2.37) 

In this subsection, we shall use the notation q = (q, q = uj) and the metric will be 
(— 1, 1, 1, 1). Although this identity cannot fully determine the vertex, it is a strong con- 
straint. In the limit q —>■ 0, 

G-\k + |) - G-\k -§)«£ %d,G~\k) + 0(q% (2.38) 

and the vertex can be chosen to be 

T»{k+ q -,k- q -)=d,G-\k) (2.39) 

In our model, the T = Green's function is given by 

G~ l {k) = g~\a)^e i ^ n{kQ \k - e k + t5k o y- a (2.40) 

The vertex would be 

q q (1 — a)e i ^ LSign ^°' ) u; a 

T ° {k + 2> k ~2 ) = g(a)(k -e k + iSk o r ^ 

r*0(fc + -£;--)= [ ^^(k) (2.42) 
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However, the sign of the imaginary part of G l (k) changes when k = 0. This leads to a 
5-fuction contribution to T° but has no influence on P^°. So, we have 



g^o v 2' 2 J g{a) 



(1 _ aW^fsigntfco) a7r 
— — — h 2i sin — - £ k 



l-a 



2' 2 ; ^(a)(fco-e k + i<Jfco)« 



u*(k) 



(2.43) 
(2.44) 



Note that the vertex functions are singular and the system is always strongly interacting. 



III. PERTURBATIONS 

In the present section, we examine the effects of perturbations. The word perturbation 
is meant in a specific sense. We augment the basic model by adding various interactions 
which we call perturbations, but the solution of the resulting models is not perturbative. We 
consider two specific perturbations. The first is the interlayer tunneling in the normal state 
and the second is the in-plane BCS pairing interaction. The first perturbation is introduced 
to understand the incoherence of single particle tunneling. The second perturbation is 
designed to bring out certain features of the in-plane pairing. 

It is convenient to use the language of functional integrals The amplitude in imaginary 
time, Z, is given by 

Z ~ J ViPlViP a e- So[r °^ ] - Sl[r ^\ (3.1) 

where So is 

Sq = ~\ ^ S^^w^CW^toW, (3-2) 

where Q(k,cu n ) was derived earlier, and the spin index is denoted by a; f3 is the inverse 
temperature, ip* and ip are the Grassmann variables representing the fermions, and the 
Fourier series 

^W = lE e ^k) ( 3 - 3 ) 

relates the imaginary time Grassman variables to Grassman variables as a function of Mat- 
subara frequencies. The various perturbations are contained in Sj. 
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A. Interlayer tunneling in the normal state 



Consider two identical layers, % = 1,2, that are coupled by a tunneling matrix element, 
t±(k), assumed to be real. Then, 



The quantity Si is given by 



SimM = - Y, t±(k) [^M^M + (i - 2) 

k,<J,U) n 

Now, So + Si can be expressed in a 2 x 2 matrix form: 

So + Si = ~ V $l„(uj n )M(k,uj n )-4) ka (uj n ). 



The column vector ipka(^n) is 



The matrix M(&, u;„) is 



M(k,u n ) = 



( g-\k,u n ) -t ± (k) ^ 



V 



(3.4) 



(3.5) 



(3.6) 



(3.7) 



(3.8) 



-t±(k) Q~ 1 (k,uj n ) 

The excitation spectrum is given by the inverse of the determinant, [detM]" 1 , analyti- 
cally continued to the real frequency axis. In particular, the singularities of the time-ordered 
Green's function are given by 

9{uj) 9{-uj) 



[detM]- 1 



+ 



e^n+i^uj) 2 ^ -t 2 ± (k) e-*™fi_(^) 2{1 ~ Q) -t±(kY 



(3.9) 



where 



n±(k,u) = [g- 1 ^)^] — (uj - e(k) ±i5). 



(3.10) 



There are branch points in the second and the fourth quadrants of the complex plane at 
uj = e(k) ± i5. The cuts are shown in Fig. (1). In addition to the cuts there are poles. Let 
us consider the case uj > 0. The poles are given by 
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FIG. 1. Complex plane for the singularities of the Green's function for interlayer tunneling in 
the normal state. 

e l7ra Q + {k,uj) 2 ^ = t 2 ± (k)e 2inn , 



(3.11) 



where n — 0, ±1, ±2, ±3, • • •. If we define 

uj - e(k) +i5 = \u — e(k)\e l6 , 
then the poles are the solutions of the following equations: 

D-(k,co) 2{1 '^ =t 2 ± (k), 



(3.12) 







+ 



2(1 -a) (1-a)' 



(3.13) 
(3.14) 



where 



D ± (k,u) = [g- l {a)uj^\^)\uj±e{k)\ 



(3.15) 



To solve these equations, let us divide the complex plane into four regions (I)- (IV) as shown 
in Fig. (2): > 6 > -f (I), f > 6 > (II), vr > 9 > f (III), and f > 6 > vr (IV). There are 
no solutions in regions (II) and (III), as it should be; otherwise, we would violate causality. 
The poles in regions (I) and (IV) are: 



uf = e(k) + t eS (k) 



it/2 tt/2 
± sin h i cos ■ 



1 — a 



1 — a 



(3.16) 
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FIG. 2. Breakup into regions (I)-(IV) to solve for the poles of Green's function for interlayer 
tunneling in the normal state. 



where 



tes(k) = g{a)t±(k) g{a 



t±(kY 



(3.17) 



In the limit a — > 0, these are the poles u = e(k) ± t±(k). As a increases, the poles spiral 
in and become degenerate at a = 1/2. Thus, at a = 1/2, no vestige of coherence is left 
|2~2fl . For a > 1/2, the poles fall off the Riemann sheet. The poles for u < are analogous. 



[G^u;)] -1 = 9{-lo) 



+ 



Analytically continued to real frequencies, the time ordered Green's function is, 

-t±{k) 

-t ± {k) e^n+ik,^) 1 -" 



( e-^n.ik,^) 1 -* 



V 



^ e iS ^Q + (k,aj) l ~ a 



-t±(k) 
e-^Q.ik,^) 1 - J 



\ 



(3.18) 



1. Spectral function 



The spectral function is 

1 



Aij(k,u) = — Siga(u))Im.[G(k,a})] 

71 



1 pi (v-e(k) 
« t ± (k) 13 



(3.19) 
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where Fij(x) is a dimensionless function. Its numerical evaluations are shown in Figs. (3) 
and (4); the explicit expressions are given in Appendix A. For the diagonal case, Fig. (3), 

0.5 




X 

FIG. 3. Dimensionless diagonal spectral function for interlayer tunneling in the normal state; 
(a = 0.25). 

there are two symmetrically situated peaks, and the spectral weight vanishes at e(k). The 
shape of the peaks cannot be described by Lorentzians, because, at high frequencies, they 
are easily seen to fall off as the power law uj a ~ l . Similarly, the off-diagonal spectral function 
is shown in Fig. (4). The, shape, once again, is distinctly non-Lorentzian, falling off as 
the power law uo 2a ~ 2 . We have also calulated the corresponding density of states. For the 
diagonal part, it is almost constant over the entire region between u c and — oj c . For the 
off-diagonal part, it is essentially zero over the same range. 



B. In-plane pairing interactions 

Consider now a two-dimensional layer in which the electrons interact through an attrac- 
tive pairing interaction. Later, we shall generalize to include bilayer systems. The interaction 
part, S/hCVv], is 

Si = -^ Vk,k'^kiM^-n( u -Wi)i/>-wi(w -w 2 )^'tM (3.20) 

" k,k',w,u;i,uJ2 
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FIG. 4. Dimensionless off-diagonal spectral function for interlayer tunneling in the normal 
state, (a = 0.25) 

For simplicity, the pairing interaction, V k k >, will be assumed to be separable, i.e. , 



Vfc, fc' — AfcAfc/ 



(3.21) 



Different pairing symmetries can be incoroprated by different choices of the functions 
Introducing an auxiliary complex field, c(u), and carrying out a Hubbard- Stratonovich trans- 
formation ||23|| , we can write 



Z ~ J VcVc*Vij* a Vij a e~^ c * {uj)ciu) - Sol ^^ ] - Sl 



The Grassmann variables can be integrated out by a linear shift, and we get 



Z ~ J x>cVc*e~^" c *^ c ^ +TrlnA ^ c *^' c( -^ , 



where the matrix A is 



n ki,k 2 — °ki,k 2 



' g 1 (ki,wi)(j w . 



fcc(wi - w 2 ) 



/ 



(3.22) 



(3.23) 



(3.24) 



The saddle point corresponds to c{uS) = 0, c*(u) = 0, for a; ^ 0. Then, denoting the 
superconducting gap, A^, by 
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we get 



A, = ^c(O), 



(3.25) 



(3.26) 



where the pair susceptibility, Xp(k), is 



1 



2 • 



(3.27) 



The excitation spectrum is determined by the singularities of the analytic continuation 
of [det A]' 1 , which is 

1 



[det Af 1 



(3.28) 



(s 2 (k) -lu 2 - i5y- a tu 2a g- 2 (a) + A 2 ' 

In an earlier publication |10| the phase was chosen incorrectly; the mistake was pointed out 
to us by Abrahams . Unfortunately, the same error appears to have propagated to later 
publications |25"1 , |2~6" . 



Clearly, there are branch points at ou — ±e(k). In addition, the spectrum consists of 
quasiparticle poles in the superconducting state. Consider e{k) > 0, and divide the complex 
plane into six regions as shown in Fig. (5). Let, 



(u - e(k) + i6y~ a = \u- e(k)\ e 



l-a i(l-a) 



and, 



(oo + e(k) -i5) l - a = \u + e(k)\ 



\l-a e i(l-a) 



n 3n 
-- < 9+ < — 
2 + 2 



3n n 
— < 9- < - 
2 2 



(3.29) 



(3.30) 



The quasiparticle poles are given by the solutions of the following equations: 



\u 2 -s\k)\ l -«ul a g-\a) = Al 



(3.31) 



and 



6+ + 9^ 



air 2mr 

1 + 1 ■ 

1 — a I — a 



(3.32) 



23 



Im (co) 



(V) 



_3ji ' %_ 
2 , 2 



(III) 




e(fc) 



(VI) 



(IV) 



(I) 



k) 



Re (co) 



2 ! 2 



(II) 



FIG. 5. Breakup of the complex plane into regions (I)-(VI) to obtain the poles of the super- 
conducting Green's function. 

where n = 0, ±1, ±2, ±3, • • • There are solutions only in regions (II) and (V). These are 
given by 



cW(A;)+A c yA;) e -*— 



(3.33) 



where 



A eS (k) = g(a)A k [g(a)^ 



(3.34) 



The quasiparticles are both damped and shifted in energy due to spectral anomaly. The 
effective gap A eS (k) collapses at a — 1. 

Analytically continued to real frequencies, the time ordered superconducting Greeen's 
function at T = is 



9{u) 



e- i ^fi_(k,cj) 1 - Q A k ^ 

A k -e^fi + (-k,-cu) 1 - a y 

A k 

A k -e^Q-i-^-u) 1 - 01 J 



+ 



V 



\ 



(3.35) 
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1. Superconducting spectral function 



The superconducting spectral function is found to be 



A s (k,u) = 



1 



e(k) 



(3.36) 



|A fc | 6V |A cfr (A;)|' |A eff (A;)|' 

where F s (x,y) is dimensionless. The explicit expression is given in Appendix B. In Fig. (6) 
we show the dimensionless spectral function for e(k) > 0; for e(k) < 0, the figure should be 
symmetrically reflected. When e(k) = 0, the spectral function has two symmetrical peaks, 
in contrast to two delta functions at ±A fc , as in BCS theory. As the energy increases, the 
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FIG. 6. Dimensionless superconducting spectral function, 
long-dashed lines are for y = 0, 1, 2, 3; (a = 0.25). 



The solid, dotted, dashed, 



weight of the peak on the positive frequency side grows and that on the negative frequency 
side decreases. The weight vanishes exactly at —e(k), but is finite everywhere else in contrast 
to BCS theory. The location of the peaks are the real parts of the quasiparticle energy. The 
energy dependence cannot be described by Lorentzians due to inherent power laws of the 
model. 
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IV. IN-PLANE GAP EQUATION 



In this section we shall solve the gap equation that was derived in the previous section. 
We want to demonstrate that a critical value of the coupling is necessary if we assume only 
a in-plane pairing mechanism, and therefore by itself it is not such a promising mechanism 
for superconductivity if the normal state is a non-Fermi liquid. We also want to show that 
the density of states at zero frequency is finite even though the mean field theory leads to 
a pair amplitude. This is because the density of states is proportional to the momentum 
integral of the diagonal spectral function and the spectrum still contains gapless excitations 
arising from the cut in the Green's function. In solving the gap equations, it did not appear 
important to use realistic one-electron band structures and the approximation of a circular 
Fermi surface is used. 



A. S-Wave 

Consider a pairing interaction that is attractive around the Fermi surface over an energy 
range ojd smaller than the cutoff u c ; for s-wave, 



-A. 



(4.1) 



Then, the T = gap equation is 

pX 



f UD Ao 



2 ' 



(4.2) 



For simplicity, we have chosen a cutoff prescription in which the integration over e is allowed 
to run between ±oo, but the remaining frequency integration is restricted between ±u)£>. 
Then, for a < 1/2, when A — > oo, the gap A is found to be 



A = . c9 -\a) (^-) 



l-a rl n 1 



l-a 
2a 



1 — a) ( na 
sm 



l-a/ V A , 



l-a 
2a 



, (4.3) 



where 
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*(§,§-«) UJ 



2a 

(4.4) 



is the critical value of the coupling necessary to obtain a non-zero superconducting gap, and 
B(x,y) is the Eulerian beta function. Note that this critical value vanishes as a — > 0, and 
one recovers the BCS result 

lim A = 2uj D e~^. (4.5) 

Consider now the equation for the transition temperature, T c : 

1 = — V - (4 6) 

Once again, we choose to introduce cutoffs on the Matsubara frequencies at ±cud, but let 
the integration over e run over ±oo. Thus, 

i=^r^E — - — r- a - (4.7) 

We find that the that the transition temperature, T c , is 

7T V2 2q -4a) \ A / ' v 7 

where the critical value of A is given by 



I -a) 

The existence of the critical value of A was first pointed out by Balatsy p7[ , although his 
discussion of the reentrant phase transition appears to be incorrect due to an incorrect 
treatment of the cutoff. In the limit a — ► 0, A c — > 0, and 

T c =^e^. (4.10) 

The 4% discrepancy in the prefactor in comparison to the BCS is result is due to the 



different choice of the cutoff. This result disagrees with those obtained by Sudb0 [25] and 



Muthukumar et al. ||26|| . Our calculation yields the BCS result in the limit a — > 0. 
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1. Density of states for s-wave 

In Fig. (7) we show the superconducting density of states for various values of the 
parameters; see Appendix B for details. The surprising feature is the gaplessness, although 
a pseudogap persists. The square-root singularity of BCS is gone, but the peak is located 
at the effective gap. The low frequency behavior of the density of states N(u) for a < 1/2 
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0.0 1.0 2.0 3.0 4.0 5.0 
CO / A 

FIG. 7. Superconducting s-wave density of states for in-plane pairing. The solid, dotted, 
dashed, long-dashed curves are for a = 0.05,0.075,0.01,0.25. 



is 



N(u) - N(0) 



a{2 1 ' a - 1) 
2 - a 



M a ) A 



2-a 



+ 



2 V V 1 uoj \2(l -a) '2(1 -a) 



(4.11) 



where 



^(0) _ 1 Kg - <>)) / ( AN ^ 
p sin(7ra/2(l - a)) \ 1 1 uj c 



(4.12) 



The frequency dependence of the first term is nonanalytic and the coeffecient is essentially 
proportional to a for small a; the cutoff and the gap dependences are weak. For, small a, 
the quadratic term competes with the first term. The background, N(0), vanishes in the 
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limit a — > 0, and one recovers the BCS result. The background also increases as a increases. 
Superconductivity is of course preserved because the saddle point of the action corresponds 
to broken gauge symmetry. 



B. D-Wave 

The zero temperature d-wave gap equation can be easily solved if we assume a separable 
kernel, where 



Vk,k> = — A cos 29 k cos 29k' 
and A > 0. Then the gap A k is given by A fc = A cos2# fc , and A , in turn, is 



A = u c g («) — 



sin 



ixa 



1 — aj \ A , 



l-q 
2a 



where 



pA c = 4a — 



■Kg -{a) 



udJ B I s - a 



i i 

2' 2 



The superconducting transition temperature can be found as before and it is 



rp _ ^£ 
7T 



2 2a - 4a I A 



i 

2 a 



The critical coupling, A c , is twice the s-wave critical temperature. 



1. Density of states for d-wave 



(4.13) 



(4.14) 



(4.15) 



(4.16) 



The superconducting density of states N(u) is shown in Fig. (8). The logarithmic peak 
at the amplitude of the gap is smeared out and there are states at zero energy. The effect 
of spectral anomaly increases as a increases. Consider now the low frequency behavior. For 
low frequencies, the integrals can be approximated and the result is 



N(lu) - N(0) 



a 



sin(^) 
, 3a' 



B 



3 — a l + a\ a(l + 2a — a 2 



2 



g(a)A 



(4.17) 



29 



1.5 



1.0 



3 



0.5 



0.0 



1 


1 


1 III 





















0.0 



1.0 



4.0 



5.0 



2.0 3.0 

co/Ao 

FIG. 8. Superconducting ci-wave density of states for in-plane pairing. The solid, dotted, 
dashed, long-dashed curves are for a = 0.05,0.075,0.01,0.25. 

Unlike s-wave, the low frequency exponent is unaffected by a. The background term is 

1 



N(0) = P 



a 



2(1 -a) 



(4.18) 



C. Collective mode 



In the superconducting state, the pairing interaction not only gives rise to the gap, but 
also renormalizes the vertex and produces collective modes; possible long range interactions 
are important as well. In this section, we consider only the pairing interaction to examine 
the effect of non-fermi-liquid correlations. 

In the random phase approximation and in the finite-temperature formalism, the vertex 



equation in the T — > limit is given by pi 
T s (k, k + q) = F n (k, k + q) - 



d d+1 k' 



V k>k ,T 3 g s {k' + q)t s {k', k' + q)g s (k')r 3 (4.19) 



(27^+! 

We shall use the notation q = (q, q = u) , where q is now the Wick rotated frequency. 
The three Pauli matrices are the conventional n, T2, and T3. The equation is of course a 
matrix equation. Near a collective mode, T s (k,k + q) goes to 00. The mode spectrum is 
then determined from the equation 
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T s (k, k + q) = - J ^^V k ^g s {k' + q)T s (k\ k! + q)Q s {k')^ (4.20) 

Consider S-wave pairing where Vk,k' = < and therefore f s (q) = T s (k, k + q). Close 
to the phase mode, the vertex is a superposition of r 3 and r 2 , but the coefficient of r 2 term 
diverges. Thus, the spectrum satisfies the following equation: 

r d d+l k' 

r 2 « / J^ym^Q s {k> + q)r 2 g s (k')r 3 (4.21) 

The right hand side of the above equation also contains a r 3 term, but it is negligible in 
deciding the spectrum. Substituting the expression for the Green's function 



G;\k) = 



1 ' g-\k) a x 
v a -g-\-k)j 



(4.22) 



we get 



1 . r d d+1 k g-\k + g)g-\-k) + a 2 
~ J (2n) d + l [\g(k + q)\- 2 + A 2 }[\g(k)\' 2 + A^y { - ' 

Subtracting the gap equation 

r d d+1 k 1 
1 ~ J (27T) d + 1 \g(k)\- 2 + A 2 ' ( ' ' 

we get 

7 (27r)^ [|£(A; + §)|- 2 + A 2 ][|£(A;-§)|- 2 + A 2 ] " l ' J 

In obtaining this equation, we have assumed particle-hole symmetry and have used the 
property g{k) = —g{k') when k Q = —k' and = —e^. 

In our model, the normal state Green's function is given by 

g~\k) = g-\a)uj«e i2 T^ k °\ik Q - e^] 1 "*. (4.26) 

Note that the imaginary part of g~ x {k) changes sign at k = 0. For small q (q > 0), we use 
the following approximations 
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G-\k+ q -)-g-\k- q -)^{ 



2^(a)- 1 sin(^f)cc£|£ k | 1 - Q for |fc | < f , 



iS^No-q-vHk)]for|A; |>f 
and linearize the spectrum close to the Fermi surface. We get 
2 ( an\ f \e\ 2 ~ 2a de 

= 4q Sin ^ — J j ^- 2(a)a; 2a| e |2-2a + A 2]2 

+ (% 2 + l^f ) (1 - «) 2 / | e 2 + w 2 h ^2 (a) ^2 + c ,2|i- a + A2] 2- ( 4 - 2 7) 

The limits of the integrals can be set to ±oo with impunity. The notation A implies average 
over the Fermi surface. Finally, the spectrum for the phase mode is given by 

fair\ „/3-2a 1 - 2a\ f q(a)A\~ A , 9 . 9 , ^ x 

2QSm (T) B (2^'2^)(^rj *A + (l-«)'(i + l«i-«l') = 0. (4.28) 

When a <C 1, to leading order in a, after the analytic continuation Q = iq , this equation 
becomes 



fiA-fi 2 + |v F -q| 2 = (4.29) 



There are two solutions 



4fi li2 = -ia\ 2 A ± V16|v F • q| 2 - a 4 vr 4 A 2 . (4.30) 
A propagating mode can exist only if 



n 2 a 2 



?f ° £ - w (431) 

where £o — vf/A is the conventional definition of the coherence length. For a fixed a, 
the mode is entirely overdamped as q — > 0. The damping of the collective phase mode is 
due to the coupling to the gapless single particle excitations of the non-Fermi liquid. As 
we have shown from calculations of the density of states, such gaplessness persists even 
in the superconducting state. In the presence of long range Coulomb interactions, charge 
fluctuations associated with the phase mode are energetically costly and the phase mode is 
irrelevant for superconductivity. 
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V. JOSEPHSON EFFECT 



It has been emphasized [10] that the Josephson pair-tunneling Hamiltonian (pair- 
tunneling action in the language of the present paper) used in the interlayer tunneling 
of high temperature superconductors is not consistent with Fermi liquid theory, but requires 
a non-Fermi liquid Green's function as discussed in the present paper. Such a Hamiltonian 
can be presently understood only in the spirit of an effective field theory; similar arguments 
were also briefly described in Ref . || . In other words, we construct a low energy Hamiltonian 
based on certain general principles, such as symmetry. This approach is necessary because 
we do not, as yet, have a precise microscopic understanding of non- Fermi liquid theory. 
Thus, we are unable to explicitly integrate out the high energy degrees freedom to arrive 
at this Hamiltonian. The difficulty is also compounded by the fact that interesting results 
arise only if, at the same time, we can suppress single particle tunneling between layers for 
which we can only give qualitative reasons. In particular, the incoherence of single particle 
tunneling between the layers, for sufficiently large a, (See also the section on spin-charge 
separation, Sec. VIII.) is the clue that such a term should not be included in the effective 
low energy theory. This appears natural to us because an incoherent process is not a sta- 
tionary process and its time evolution is not an oscillating exponential. The essential reason 
as to why single particle tunneling is incoherent is that the electronic quasipartcles are not 
well defined in a non-Fermi liquid, and the amplitude has to be reconstructed for an elec- 
tron to tunnel between the layers. Such a process necessarily involves overlap of iV-particle 
wave functions. In contrast, the superconducting wave function is an antisymmetrized wave 
function of Cooper pairs, and the tunneling matrix element only involves the overlap of a 
pair. It is the nature of a Bose condensed state that all pairs behave the same way. Thus, 
tunneling of Cooper pairs cannot be plagued with orthogonality catastrophe of iV-particle 
wave functions. The aim of the present section is to give reasons as to why the pair tunneling 
Hamiltonian is reasonable from the point of view of effective field theory. 

It is essential to realize that conventional Josephson effect in a Fermi liquid cannot be 
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recast in the form of an effective Hamiltonian [JTOj] . The well known expression for the critical 
current for two conventional superconductors with equal gaps |f28f, 



where i?N is the normal state tunneling resistance, is unusual in two respects. As T — > 0, I c 
is a nonanalytic function of the gap (oc |A|), and is analytic only close to T c (Note that, close 
to T c , the magnitude of the Josephson effect, derived from Ginzburg-Landau theory p9| , is 



proportional to A(A/T) and the analytic result is due to the thermal smearing.). It is also 
remarkable that I c is independent of the cutoff, uip, the Debye energy. The dependence on 
the Fermi energy is buried in the definition of R^. This is a manifestation of Fermi liquid 
theory. Consider the expression for the critical current at T = 0, given by 

/ c ocE|T k , q | 2 ^^— ^r, (5-2) 

kj q Ac + ^q 

where = ^4:(k) 2 + | Ak| 2 . The integration over the momenta can be converted to energy 
integrations which extend only to the Debye energy. If the energy integrations can be 
extended to infinity with impunity, it can be seen by scaling the energy variables that 
I c oc |A|; the correction is of order (^)- This result hinges on the Fermi-liquid scaling of 
the single particle spectral function: A(Ak, Auj) = A~ l A(k, uj). Any non zero value of the 
scaling exponent, a, will result in a cutoff dependence of the critical current J c , and if this 
cutoff is much larger than the gap, the critical current can be an analytic function of the 
gap, proportional to A(A/u c ). This allows us to reformulate Josephson effect in terms of 
an instantaneous pair tunneling hamiltonian on energy scales smaller than the cutoff. 



In fact, explicit calculations JTOj using the non- Fermi liquid spectral function of the 
present paper leads to 

/ c cx(|T kiq |> F 2 A(-)/(a,-), (5.3) 



where the angular brackets imply an energy average and f(a, — ) is a complicated function. 
When £■ « 1, /(0, A/u c ) = u c /A, and we recover the Fermi liquid result. For a > 1/4, the 
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dependence of /(a, — ) on disappears entirely in the limit ^ < 1, and it is a function of 
a alone. A closed form expression in the regime ~ > a > ■? was noted earlier |1(J. For a not 
too small, but less than |, the dependence on is weak. 

We now construct a pair tunneling Hamiltonian. This Hamiltonian is certainly justified 
for a > j, which is also the threshold at which incoherence sets in in the spin-chrage 
separated model (See Sec. VIII); for the case without spin-charge separation the threshold 
of incoherence is a — \. In the tunneling between the layers, the transverse wave vector of 
an electron is conserved, and this conservation is an exact symmetry. Thus, consider the 
effective action 

Sj = J dt'dtd^d^gQ^, k 2 , |t - t'\) [tfiWVtfjS.CO^W^CO + (1 - 2)] . (5.4) 

In view of the discussion above, we may set g(k 1; k 2 , \t — f\) = 9 ^' k2 ^ (£ — t'). Because in a 
non- Fermi lquid the tunneling process is so widely and uniformly spread over frequency that 
the Kernel is instantaneous. If we now repeat the scaling arguments of Sec. HE, we see that 
the scaling dimension of this term at the non- Fermi liquid fixed point is s _< - 1_2a ^. This term 
is therefore relevant, diverging as s — > 0, as long as a < \ and becomes marginal at a = |. 
In contrast, we have seen in Sec. HE that the conventional BCS pairing interaction (the 
in-plane mechanism in the interlayer tunneling theory), which involves integration over three 
independent wave vectors, is irrelevant, scaling as s 2a , hence vanishing as s — > 0. Indeed, 
we have seen in Sec. Ill that a coupling larger than a critical coupling is necessary to drive 
a superconducting transition, if we use a BCS pairing kernel. 
It is reasonable to assume a separable form for g(ki, k 2 ): 

^(k 1 ,k 2 ) = -t ± (k 1 )t J .(k 2 ) (5.5) 

where £j_(k) is the single electron tunneling matrix element that can be determined from 
microscopic band structure calculations. An expression for ij_(k) based on general consid- 
erations of symmetry and analyticity was proposed in Ref 0; it is 

t±(k) = (cos k x a — cos k y a) 2 . (5.6) 
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This expression for ij_(k) is now confirmed from more detailed considerations of actual band 
structure computations p0| . 

As in the original BCS theory, the pair tunneling action can be simplified, and a reduced 
action can be used for computations within the mean field approximation. This reduced 
Hamiltonian, which involves only pairing of states (k f) with (— k J,), is 



Note that it is misleading to consider the relevancy of this reduced action with the dimension 
of the fermion operators determined from the fixed action of the normal state. The reduced 
action is equal to the original one only in the BCS subspace. 

In addition to pair tunneling, we may also construct a particle-hole pair tunneling term 
given by 



sr h = J dmdVk^, k 2 , it - 1>\) [^w<Uo<UovCJUo + (i - 2) 



Similar to the single particle tunneling, the importance of this term is difficult to decide 
because this process is likely to be highly incoherent in a non-Fermi liquid for precisely the 
same reasons. In other words, the time dependence of the kernel may not be replaced by 
8(t — t'). If, on the other hand, the system has a tendency to form a spin density wave, this 
term may be more important than the pair tunneling term. 



For interlayer tunneling in which there are both in-plane attractive interactions and 
Josephson pair hopping between the planes, we take So to be 




(5.7) 



(5.8) 



VI. INTERLAYER PAIR TUNNELING 




(6.1) 



Similarly, 
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s)- p 

and 



- r 



k,k' 



~ UxWX^iu - u 2 )^>(uj 2 ) + (1 - 2) (6.2) 



X 1 ) 



(1), 



-■gs E W^fC^)^!^- a; 1 )^(^-^)Vk 2 T ) (^) + (l-2)] (6.3) 



The term S 1 } -13 corresponds to a BCS in-plane pairing kernel as discussed in Sees. Ill and 
IV, while the term Sj is the interlayer pair tunneling term introduced in Sec. V. This model 
was discussed in some detail in Ref. [[|, where, for simplicity, the single particle spectrum 
was assumed to be that of a Fermi liquid. The assumption was that for in-plane motions the 
quasiparticles are recovered in the superconducting state. This assumption is not entirely 
justified. As we have seen in the previous sections, the spectral function is in fact quite 
broad; see Fig. (6). Thus, in the present section we include properly the non- Fermi liquid 
features. 

The Hubbard- Stratonovich transformation for the combined action is not entirely 
straightforward due to the local- in-k pair hopping term. First, let us rewrite the Josephson 
term as 



S J j = 2 [MK u)*A(k, uS) - R 2 {k, u) + J 2 (k, 



(6.4) 



k,ui 



where we have ignored an unimportant single particle term that can be incorporated in the 
band structure. Here, 



* 1 /Tj(k) 

=2\/^ 



E[V'k 1 T ) *M^>-^i) + (i-2; 



(6.5) 



I rak) 

4V P 3 



tf( k ><") = 7^/^ E KV* M^ll^-u,) -(1 -> 2)] (6.6) 



Ak,^) = - 



1 Zj(k) 



4 V f3 3 



E [4f - u,) - ^l(o, - u^'M - (1 



',(1)/ 



(6.7) 
(6.8) 



The in-plane terms are decoupled, as usual, by using 
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4'M = ^E A(k)4 1 T ) *(^)^(c - Wl ) 
4» = i E MV^> - ^ 

P k,u> 

We also include in the functional integral the factor 



(6.9) 
(6.10) 



I]cxp [-Z^M^iM - B*{u)B 2 {u)\ , 



(6.11) 



and the factor 



J|exp [-2A*(k,w)A(k,w) -2R 2 (k,u) -2I 2 (k,u) 



(6.12) 



and integrate over f?i(u;), £*M> B 2 ((v), B%(u), A(k,uj), A*(k, u), and over two real fields 
R(k,uj), and 7(k, u). A set of linear shifts 



s 2 M -> s 2 m - 4M 

A fe -> A fc -i(k,w) 
.Rfc — > -Rfe — ii2(k, a;) 
4 -> 4 - /(k, w) 



(6.13) 
(6.14) 
(6.15) 
(6.16) 
(6.17) 



eliminate all the quartic terms. The fermions can now be integrated out and the action can 
be written as 

s = E fo'M^iM + 5 2 *M£ 2 M] + E [^*( k > ^M( k > u) + R 2 (k, u) + / 2 (k, w )' 

-InDet [M (1) M (2) ] , (6.18) 

where and M^ 2 ^ are two 2x2 matrices in k-uj space. The problem can be simplified if 
we note that at the saddle point 

A(k)2M«) + f ^(k, u) + ^ + ^ = 



\(k)B* 1 (u) + JT J (k)[A*(k,u;) + 



iR(k,u) I(k,cu) 



(6.19) 
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Therefore -R(k, uj) = I(k, to) = 0. At the saddle point the determinant of the product of two 
matrices simplifies to 



Det 



M (1) M (2) 



1 



^- 1 (k,^)e;- 1 (-k,^-^ 1 )+/3- 1 |A(k) J B 1 + v /T J (k)A(k)| 



^- 1 (k,^)e;- 1 (-k,^-^ 1 )+/?- 1 |A(k) J B 2 + v /T J (k)A(k)| 



, (6.20) 



where all the fields are evaluated at u = 0: Bi = Bi(0), B 2 = -82(0), and A(k) = A(k, 0). 
On taking functional derivatives, we get, 

SS 







SS 



1 



8B{ 



A(k) A(k)Sx + JTj(k)A(k) 



P k*, g-^K^g-^-k, -00) + p- i \\(k)B 1 + JTj(k)A(k)\ 2 



8A(ky 



2A(k)--Y: 



Tj(k) A(k)5 1 + JTj(k)A(k) 



P g- l (k,u)g-\-k, -u) + P-i\\Ql)B x + ^(^^(k)! 2 



A(k) 


\{k)B 2 + y/Tj(k)A(k) 




g-\k,u)g-\-k 


-u)+f3-i\\(k)B l + v 


'Tj(k)AQ£)\* 



n - SS - B 1 V 



The saddle point is at B\ = B 2 = B. Now if we define the total gap A k by 



(6.21) 



(6.22) 



(6.23) 



A k = X(k)B + JTj(k)A(k) 



(6.24) 



then 



A k [1 - x(k)T,(k)] = A(k) £ A(k')A k , X (k') 



(6.25) 



The pair susceptibility in terms of the total gap is the same as defined earlier. Note that the 
only difference from the result in Ref. || is the presence of non-Fermi liquid effects in the 
pair susceptibility. The derivation may appear to be a little more involved than necessary, 
but the formalism that we have set up could be quite useful if one is interested in the 
theory of collective modes |3l| . This equation can be solved for any symmetry of the pairing 
interactions. 

At T = 0, this gap equation can be easily solved on the Fermi surface. Consider the 
anisotropic s-wave case 0: A(k) = VX At the Fermi surface 



A kF = A + Tj(k F )A kF x(kF) 
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(6.26) 



where 



A = A$> k , x (k') 



(6.27) 



First consider the case a <C 1/2. Then, we can let the cutoff to oo, because the integral 
converges. We get, if the Josephson pair hopping term dominates, 



A 



u c g l {a) 



~Tj(k F y 


l-a 
l-2a 


1 


1- 
1- 






[2(1 a) sin 2( £ a) _ 





(6.28) 



This, as noted earlier | ]25(| , is encouraging, because local-in-k interlayer hopping produce a 
gap on the Fermi surface and there is no need for a critical coupling, necessary for purely 
in-plane mechanism discussed earlier. As a — > 1/2, the above approximation of neglecting 
the cutoff uid cannot be used. However, the inegral can be easily approximated with the 
cutoff in place. The result, for a — > 1/2, is 



A kF S3 



1 / \ I U D 

u c g {a) ( — 



1 7TLJc(l — as) 

P. Tj(k F ) 



(6.29) 



Note that the precise results are not in agreement with those of Ref. ||25|| . This is due to an 
incorrect choice of the phase of the Green's function. 

The density of states as to —>■ can also be estimated, and surprisingly the behavior 
found earlier for the in-plane pairing is unchanged. We find for a C 1/2 and uo — ► 0: 



N{uo) - iV(0) _ /u\ 2 - 2a / 2 -2a 
p \xj V 2 — a 



(6.30) 



where A < x < A max ; A max = max[A c fj(k)], and A = min[A e fj(k)]. Thus, there is a finite 
density of states at uj = 0, and the low frequency behavior is u 2 ~ a . 

It is interesting to compute tunneling density of states with a realistic Fermi surface 
corresponding to a two-layer high temperature superconductor, such as Bi2212 or YBCO, 
assuming an s-wave in-plane pairing kernel. For this we choose |32 



Ek = — 0.5(cos k x a + cos k y a) + 0.45 cos k x a cos k y a, 

where a is the lattice constant. The pair-hopping term Tj(k) was chosen to be 
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(6.31) 



T 

Tj(k) = — (cos La — cos k v a) . 
16 



(6.32) 



The other parameters are as follows: Tj = 0.05 eV, Ef = —0.315 eV, Ao = 0.01 eV, and 
io c = 0.15 eV. The results are shown in Fig. (9). The tunneling density of states are 
qualitatively similar to those of Renner and Fischer |33| . 
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FIG. 9. Superconducting density of states for interlayer tunneling with realistic band parame- 



ters. The solid, dotted, dashed, long-dashed curves are for a = 0.05,0.075,0.01,0.15. 



VII. SPIN-CHARGE SEPARATION 



In this section we shall consider a simple model for spin-charge separation [33|. The 
time-ordered Green's function at T = is 

G(k,u) = — , (7.1) 

lo^(uj — v s k + i6up (u> — v p k + iSup a 

where v s and v p are the spin and charge velocities (v p > v s ), k is measured from the Fermi 
surface, g(a) is the same as before, and uj c is the cutoff. The corresponding spectral weight 
function is 

A(k,tu) = -\lmG(k,u;)\ (7.2) 

7T 



41 



(7.3) 



^\G(k,u)\sm(™), when (u - v p k)(cu - v s k) > 
^\G(k,uj)\ cos(^), when {uj - v p k){uj - v s k) < 

It has two power-law singularities and high frequency tails. The spectral weight in the region 

between the two singularities is larger by a factor of cot( z y). 

The momentum distribution function can be obtained by the formula n(k) = 

J° UJc A(k,uj)duj and yields the value 1/2 for n(0). At finite k, and for v p \k\ « uj c , we 

get 



n(k)^ 1 --^ 



[v p - v s )k 



v.. 



OJr. 



(7.4) 



where the function f(a,y) is given by the integral 

dx 



f(a,y) = a [ — 

JO J2 



(x + y — 1) 2 

We recover the earlier result for n(k) without spin-charge separation when we set v p — v s , 
which results in f(a, 1) = I. The density of states, n(u), related to the spectral function by 
the formula n{uj) = J2k w), i s constant when v p and v s are independent of uj and k. 

If we consider, once again, the problem of tunneling between two planes, we find that 
the inverse of the matrix Green's function, G t (k,cu), is 



[G t {k,uo)}- 1 



( G~\k,u) -t ± (k) 



(7.6) 



^ -t ± (k) G-\k,u) 

The quasiparticle poles are given by the equation det \Gt l (k,oj)\ = 0, more explicitly by 



CO 



2 c a g- 2 (a)e ta7TSien{uJ \uj - v s k + iu5)(oj - v p k + iuj5) l - 2a = t\(k). 



(7.7) 



From the phase and the amplitude of this equation, we get 

Sign(uj)an + 6 S + 9 p (l - 2a) = 2Ntt, 

uj 2a \uj - v s k\\^ - Vpkl 1 ' 20 = g 2 {a)t\{k). 



(7.8) 
(7.9) 



For k > the equation for the phase restricts the poles in the shaded area shown in Fig(10). 
In general, there are two poles ui^ and , where Re u£ > v p k, Im oj£ < 0, and Re oj^ < v s k, 
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FIG. 10. Complex plane depiction of the singularities of the Green's function for the spin-charge 
separated model. 

Wfelm ojJ, < . However, when a > \ and t±(k) is less than a critical value, the pole Uk+ 
is absent, which clearly indicates supression of coherent tunneling. This is an important 
respect in which the spin-charge separated case differs from the previously discussed case of 
no spin-charge separation. To understand why the pole disappears, consider the equation 
for the phase: 

cox + 9 S + e p {\ - 2a) = 0. (7. 10) 

Because 9 p > — |, 9 S is bounded by 9 S < |(1 — 4a). When a > |, this inequal- 
ity confines the pole in the shaded area as shown in Fig(ll). For oo belonging to the 
shaded area, the left hand side of the amplitude equation has the minimum: u^ a (v p k — 
v s k) 2{ - l - a) \ sec(0(a))||tan(0(a))| 1 - 2a , where <p(a) = f(l-4a). For g 2 (a)t\(k) less than this 
minimum value, the equation for the amplitude can not be satisfied. So the pole Uk+ does 
not exist for a > \ and t±(k) less than a critical value given by the minimum. 

We now discuss the superconducting properties. The calculations are similar to those 
discussed earlier. For S-wave pairing, the gap is 

A = g~ 1 {a)u D 
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L — (Jt 



1 - 



A 



(7.11) 
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FIG. 11. The restriction of the pole in the shaded region of the complex plane. 



where 



an 



pg 2 (a)h(a,^) \u D J 



2a 



h(a, f3) 
h'(a,(3) = 



sec 2a 9d9 



cos 2 6 + f3 2 sin 2 6 
d6 



>o [cos 2 9 + (3 2 sin 2 0]s=5= 
For rf-wave, the gap is given by = A cos 9k, where the amplitude A is 



A = g 1 (a)u D 



2a 



1-a 



(7.12) 



The critical coupling A' c = 2A C . The equation for the transition temperature, T c , is unchanged 
from the case without spin-charge separation, provided the present definitions of the critical 
couplings, A c and A' c are noted. 

We find that the quasiparticles are damped, and the real part of their energy is greater 
than \v p k\. Close to u = 0, the density of states, N(u), is smaller than that of the normal 
state. For s-wave pairing, we get 

2-a 2-3a\fg(a)A [vp~\~ 



N(0) = p 



1 - (l-a)B 



2 -2a 2 -2a 



(7.13) 
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and 



N(u)-N(0) \u 



2-2a, ,a 



P 



oc a- 



For (i-wave pairing, we get 



N(0) = p 



1 - -B 

2 



1 2 — 3a\ ( g(a)A fv~ p 



2a 2 - 2a 



and 



N(u)-N(0) \u\ 
oc 



(7.14) 



(7.15) 



(7.16) 



P A 

As before, the linearity of the ci-wave density of states is unaffected by a. 

We conclude that the spin-charge separated model behaves qualitatively similarly to the 
model without spin-charge separation. However, we found one important difference between 
the two models with respect to the coherence of tunneling between two layers in the normal 
state. In the spin-charge separated model, the coherence is absent if a > 4, and the tunneling 
matrix element t±(k) is less than a crtical value. Given the anisotropy of t±(k) @, this can 
lead to interesting results. 



VIII. CONCLUSIONS 

In the present section we shall summarize our work and point out the directions in which 
our theory must be improved and extended. We have shown that the notion of spectral 
anomaly is a powerful one. It allows us to formulate a framework with which we can un- 
derstand non-Fermi liquid effects on superconductivity, regardless of the microscopic origin 
of the spectral anomaly. We believe that any low energy effective field theory of gapless 
fermions can be characterized by an anomalous spectral function that reflects the critical 
state of the fermions. The results obtained from the spin-charge separated model are quali- 
tatively similar to those obtained from the model wthout spin-charge separation. However, 
the single particle tunneling between the layers shows a greater degree of incoherence in 
the spin-charge separated model. We have also shown that the concept of spectral anomaly 
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provides a basis for understanding the features of the interlayer tunneling theory of high 
temperature superconductors. 

Spectral anomaly has two important effects. It makes the the single particle spectrum 
incoherent while retaining some of the features of a Fermi liquid. The spectrum is gapless 
and therefore the system responds as a metal: an infinitesimal electric field generates a 
conducting response. We have also seen that one can still define a surface in the momentum 
space at which the nature of the excitations changes, as at the Fermi surface of a Fermi liquid. 
However, the momentum distribution is not discontinuous at this surface but its derivative 
is singular. Spectral anomaly has also the negative effect that conventional BCS pairing 
is irrelevant (in the renormalization group sense). We have seen that one needs a coupling 
larger than a critical coupling to achieve superconductivity. Nonetheless, the interlayer pair- 
tunneling Hamiltonian is relevant (not marginal). This is particularly interesting because 
many of the features are quite distinctive and reflect the dominating effect of the interlayer 
tunneling rather than the small residual in-plane attractive mechanism, which can of course 
have any symmetry. According to the present analysis, the c-axis tunneling density of states 
is finite at zero frequency, but a well-defined peak still persists at a finite frequency unless 
the anomaly exponent, a, is very large. While such effects are reported, it is not clear if 
these are experimental artifacts or not. However, if such effects are not artifacts and can be 
measured precisely, we can infer the value of a from such experiments. 

In the present paper we have been unable to discuss the renormalization of the composite 
operators such as current. Although we have used the Ward identity to indicate how the 
vertex should scale, it has not been possible to go further. In our opinion, it is not possible to 
do so without more powerful methods, such as the unitarity expansion formulated many years 
ago by Gribov, Migdal, and Polyakov The renormalization of composite operators is 
important because most experimental probes couple to composite operators. In this respect, 
it is necessary to go beyond the leading scaling forms, as the dependences on temperature 
and frequency are often determined by the the presence of dangerously irrelevant operators 
as argued by Nayak and Wilczek |36| . 
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Our justification of the interlayer tunneling Hamiltonian consisted of a two step process. 
We first showed that the single particle tunneling is incoherent in a non-Fermi liquid for 
sufficiently large a, which need not, however, be much larger than | to 4. We then argued 
that because of incoherence such a term should not be included in the effective low energy 
theory. Renormalization group arguments, similar to those discussed in our paper show, 
however, that such a term is relevant as long as a < 1. Josephson pair tunneling was then 
argued not to be severely affected by the spectral anomaly and can be cast into an effective 
low energy Hamiltonian. It can be argued that it should be possible to begin with a single 
particle tunneling Hamiltonian and arrive at a low energy theory in which single particle 
tunneling is absent but pair tunneling is present as a result of correlation effects that give 
rise to spectral anomaly. Our inability to do so is intimately intertwined with the neglect of 



vertex corrections in the tunneling process as was remarked by Wen [11 



Ultimately, it will be necessary to provide a microscopic derivation of the spectral 
anomaly. Beginning with the pioneering work of Anderson many authors have inves- 
tigated the notion of tomographic Luttinger liquid in two dimensions. Similarly, many au- 
thors have focused on an underlying gauge theory of strongly correlated electrons that could 
possibly give rise to non-Fermi liquid behavior. Recently Khlebnikov [I2~| has approached 
interlayer pair tunneling from this perspective. Such ideas seem promising to us. 
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APPENDIX A: SPECTRAL FUNCTION I 



We derive here the spectral function corresponding to one particle green's function of 
the interlayer tunneling problem discussed in Sec. ( |III A|) . 

Aij(k, uj) = -isign(o;)Im [G(k, uj% , (Al) 

where for both the regions e(k) < uj < and uj < e(k), uj < 0, 

sin f Sf] D~(k, uj) l ~ a \t\(k) + D~(k, uj) 2 ^' 

lm[G(k,uj)l, = ^ L 2 — — l —. (A2) 

L 11 [D~(k, u;) 2 ( 1 ~ a ) — cos(aiT)t 2 ± (k)] + sin (a7r)t\(k) 

For the regions uj > 0, to > e{k) and < to < e(k), Im [G(k, u))] n is negative of the expression 
given in Eq. ( |A2| ). In the more compact form, the spectral function is 

^• w) =d^ F »(w)" (A3) 

where the dimensionless function Fu(x) is given by 



sin (f ) 


\x 




X 


2(1-")] 


7r[|a;|4(i-«) - 


- 2 cos(a7r) 


I 


2(l-a) + 1] 



Fl1 ^ vrflxl 4 ^) -2cos(avr)|x| 2 ( 1 ~«) + 1] ' ^ 
Similarly, the spectral function related to tunneling between the layers is related to the 
off-diagonal element by 

A 12 (k, uj) = --Sign(cu)Im [G(k, w)] 12 , (A5) 

where, for e(k) < uj < 0, 

T m sm{an)D-{k,uj) 2 ^- a k ± (k) 

Im \G(k, uj) , 9 = — ^ — . (A6) 

L Jl2 [^-(A;,^)^ 1 -") -cos(a7r)t 2 L (A:)] 2 + sin 2 (a7r)tl(fc) 

In the region uj < 0, uj < e(k), Im [Gr(fc, o>)] 12 is negative of that given in Eq. ( |A6|) . In the 
region uj > 0, uj > e(k), Im [G(k, uj)] 12 is negative of that given in Eq. ( [Alf ). Finally, in the 
region e(k) > uj > 0, Im uj)] 12 is the same as that given in Eq. ( |A6| ). 
Therefore we have 

A ^m F -(^)- (A7) 

where 

Sign(x) sin(a7r)|x| 2(1 - a) 
121 J 7T [[Irrl 2 ^) - cos(aTr)] 2 + 1] ' { } 
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APPENDIX B: SPECTRAL FUNCTION II 



It is interesting to calculate the superconducting density of states and the spectral func- 
tion. First, consider the spectral function. The imaginary part of the matrix element, 



[G?{k, 



UJ 
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, of the retarded Green's function is related to the spectral function by 



A 9 (k,u) = Im Gf(k,u 



IT 
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In terms of the time-ordered Green's function, 



A s (k,u) = Sign(w)Im [G s (k,u)] u 

7T 



(B2) 



First, consider first uj < 0. There are three separate case: (1) e(k) < uj, (2) ui < e(k) < —uj, 
and (3) — uj < e(k). 

1. e(k) < uj 

In this region, we have e(k) < uj < < — uj. Therefore, both uj — e(k) > and 
-uj - e(jfe) > 0. Then, 



Im [G s (k,uj)} n 



sin 



D+(k,uj) 



l-a 



[D+{k,u)D-{k,uj)Y 1 -^ + Al 



(B3) 



2. uj < e(k) < -uj 

In this region, there are two possibilities, either uj < < e(k) < —uj or uj < e(k) < 
< — uj. If we remember that we go counterclockwise around one branch point and 
clockwise around the other, we get, for both of these cases, 

sin (eg) D + (k, uj) 1 -* \A\ + [D+(k, u)D~(k, cu)}^' 



Im [G 8 (k,u)] n = 



[[D+(k, u)D-(k, lo)]^ - cos(a7r)Al] 2 + sin 2 (a7r)A 



(B4) 



3. — uj < e(k) 

In this region, both uj — e(k) < and —uj — e(k) < 0. Therefore 

sin(^) D+(k,ujy- a 



Im [G s (k,u)] u = 



[D+{k,uj)D-(k,u)]( 1 - a ) + A 2 k 



(B5) 
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The case u > can be treated identically and we give the results below, written in a 
more concise form. 



1. \e(k)\ > to. 



Im [G s (k,u)] 



sinf^f )D + ( y k,u 



,1-a 



11 



[D+{k,uj)D-{k,uj)]^-<*) + A 2 



(B6) 



2. |e(fc)| < w. 



Im [G s (k,u)] n 



sin (<f ) £>+(£;, a;) 1 "" [A|, + [D + (k, uj)D~(k, cj)] (1 ~ q ) 



[[D+(fc,w)D-(A;, c^p"*) - cos(a7r)A 2 ] 2 + sin 2 («7r)A| 



(B7) 



In the scaled form, the spectral weight can be written as 



A s (k,u) 



1 



e(k) 



(B8) 



where 



sin (f ) 



7T 



f (|y| - 


x\) 


\x + y\ 1 ~ a t 0(\x 


- | 2 /|)|x + 2 /| 1 - a (l + 


x 2 - y 2 


l—a\ 




x 2 — 


y 2 


+ 1 [|a; 2 - y 2 


1 a — cos(o;7r)] 2 + sin 2 (o;7r) 



(B9) 



The BCS result, for comparison, is 



A BCS (k, u) = u 2 k 5(u - E k ) + v 2 k 8(u + E k ), (BIO) 
where E k = \Je 2 (k) + A 2 ., and 



2 1 /\ , 



M -2 1 + t ' (BU) 



.2 _ 1 A 



^=2 rir • (m2) 



The superconducting density of states, N(u), is 



iV(cu) = --Sign(cu) ]T Im [G S (A;, u)] u . (B13) 
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